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1.0 INTRODUCTION 


1.1 Alms of the present work 

The aim of the work described in this recart is to adapt or develop 
a jet- i n-cro5sf low itmdel for estimating wind tunnel interference with 
primary emphasis upon non- Imp i ng i ng cases. A cor respond i nq computer 
program Is required which can be used on an optional basis with the Part I 
program for tests Involving mode 1 with llftinq jets O’* fans. 


1 . 2 Bac kqround 

Part I of this report describes a method for calculating wind tunnel 
blockage and anq le-of-at t ack corrections from pressure measurements made 
along the tunnel surfaces. Theoretical flow models employed are consiructed 
from line singularities which represent the model and it's wake as 'perceived' 
by the tunnel. The flow models can accommodate wing sweep, anq le-of-at tack 
and offset from the tunnel center, though in many cases it has been found 
that an unswept planar model suffices. The approach may be considered semi- 
inverse Since the details mentioned are user- spec i f i ed yet the singularity 
strengths and their axial locations are determined from wall pressures using 
influence matrices. 

It would be very useful If the above approach could be extended along 
the same lines to lifting Jets and a 1 K>w traiectory shape, for example, to 
be implied from the wall pressures. Hcsvever, on reviewing the logistics 
needed for doing this - particularly the likely size and conditioning of 
the influence matrices concerned - It soon became apparent that a more 
explicit approach was preferable. Subsequent experience, described later, 
has shown that use of an explicit model for the jet plume is probablv 
essen t i a I . 

It appeared Initially that, with the exception of wall pressure 
measurements, most of the needed experimental data and flow modeling exper- 
ience would be available in existing literature. However an extensive 
review of experimental data revealed that most experiments concerned je’ts 
emergent from a plane or from models which are relatively complex. Jel- 
from-pipe data, which was considered more suitable, was relatively scarce. 

New experiments were therefore p 1 anned , us i ng jets-f rom-pipes, to determine 
wall pressures and to find parametric relationships between the new jet- 
from-pipe and the existing je t - f rom-p 1 ane flow field data. 

A review was made of the very many theoretical models which have been 
proposed for the jet- i n-crossf low. Most had to be rejected because thev 
are too complex for the present application. However the models proposed 
by Fearn (Ref. 1 ) . by Heyson (Ref. 2) and by Williams and Wood (Ref. 3) 
Included at least some of the necessary physics and were reasonably simple. 
These were therefore prime candidates in the earlier stage of the work and 
it was hoped that one or mc>re could be used directly. 
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The principal t U'lw model property required for the present work is 
the ability to predict far-fteld flows. A good test of this is to attempt 
to predict tunnel wall pressures. It will be seen later that all of the 
candidate methods failed this test quite badly. In retrospect, it is 
apparent that all were suitable only for near- field app 1 i cat ions . 

As originally perceived, the main tasks were to interpret and 
organize existing experimental data and modeling techniques for tumiel inti'r 
ference prediction. It has turned out that significant new studies were 
needed in both areas In order to build upon existing technology. These 
studies are the subject of the rest of rhis repv'»rt . 


1.3 Layout of the Present Report 

Section 2 will describe, in broad term^, the way in which jet-in- 
crossflcxv and "rest of model" interference effects are determined and 
combined. Sections 3 and h concern test details and results for jet-from- 
p i pe tests conducted as part of the present study. Jet hardware and 
calibration procedures are described in more- t han- usua 1 detail because thf 
extra attention given to these aspects paid off well. The development of 
a new flc>w model, designated the 'VSD' model for Vortex/Source/Doublet, is 
described in Section S and comparisons are made with predictions using 
other methods. Section 6 gives the fairly comprehensive conclusions 
which arose from the present work. 

Ill the interest of shortening the main part of this report, the 
majority of the experimental traverse data Is presented as Appendices A and 
0. For similar reasons, comments concerning Hevson's method for interfer- 
ence estimation appear in Appendix C. Program documeiUat ion and listings 
are given in Appendix D. 



2.0 


interface to the part I PROGRAM 


For tne reasons mentioned above, the jet- i n-crossMow modeling program 
developed here is used as a pre-processor to the main. Part I program. 

Figure 2.1 depicts the general principles of operation of the Part I program. 
Prior to a test, a theoretical flow model is constructed using model span, 
sweep and ang 1 e-of-at tack details, as needed. An influence matrix derived 
from this model is used in conjunction with wall, roof and floor pressure 
data to determine source and vortex strength distributions as a function of X. 
Tunnel-induced ang 1 e-of-at tack and blockage increments are then calculated. 

For jet- or fan-lift models the general approach is the same hut the 
above routine is preceded bv explicit j et- i n-crossf low calculations (see 
Figure 2.2). These accomplish two objectives. Firstly, a direct estimate 
is made of the distributions of Au and A.'i, for subsequent addition to 
corresponding • rest-of-mode 1' data. Secondly, the effect of the jet-in- 
crossflow (or jets) is removed from the measured wall pressure signatures 
prior to further processing. This step is very important because it removes 
substantial non-planar effects from the data which would otherwise cause 
errors due to 'cross' effects. In particular, offset blockage effects of 
a strongly penetrating jet would be returned eventually as tunne 1 - • nduced 
angle-of-attack in the absence of the preprocessing step. 
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3.0 test rigs and procedures 


3.1 Introduction 

Figures 3 1 and 3.2 give a general view and dimensions of the test rig 
for the jet-in-crossflow experiments. The main test measurements comprised 
in-plume three-component velocity traverses, using the rake of 5-holed pitch 
yaw probes shown, and wall pressures along the test section walls and roof 
(see subsection 3-^). Boundary layer tangential blowing was available at the 
roof, from a location 13 -lnches ahead of the jet exit, and was applied when 
surface pressures on the roof centerline indicated the presence of flow break 
down. Sufficient blowing was applied to remove the standing vortex, ahead o 
the jet impingement point, when it occurred. This parallels the approach 
discussed in detail in Part I of the present report for jet-flapped wings. 

In practice, the tangential blowing was needed only with the 3”'nch jets at 
high jet velocity ratios. 

Figure 3.3 tabulates the jet configurations used and typical mainstream 
and jet velocities for the nominal test velocity ratios. With the 3-lnch 
jets at velocity-ratio 8, a reduced mainstream speed was necessary because 
of supply limitations. 


3.2 On-line Tunnel Blockage Corrections 


The conventional solid-plus-wake blockage equation has the form 

Uc - Uo (1 + ^w^ 


where U^. is the corrected velocity at the model location 

U is the calibrated empty tunnel velocity at the model 
^ 5 * ‘w solid and wake blockage factors. 


The above form was used because the full wall pressure signature approach 
cannot be implemented without the flow model which will be developed from 
the present experiments. A predecessor of the pressure signature method, 
formerly designated the 'q-pot' approach, was used to determine 1^- ^ con 

ventional calculation was used for Cj, the solid blockage based upon jet 

dimensions. The velocity increase AU , due to wake blockage. i> infer ed 


pipe 
V i a 


via linear assumptions from the pressure decrease .^Cpt , between the con- 
traction exit piezometer ring and the (atmospheric) breather slots at the 
start of the first diffuser. Thus 


c 


w 



1 

2 


AU 

U 


X 


m oo 




(3.2) 


On supplying the necessary calibration constants and removing the empty- 
section value of ACp . a working equation for blockage-corrected velocity is 
obtained at the model location as 
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(3.3) 


. i - 1.0182 (I + Cj + 0.25 (.ICpt - .04903)1 4p 

- 1.0182 (q^ / q^) ip 

where ip Is the measured contraction pressure drop 
and ACp^ « / 1 .0182 Ap. 

The corresponding corrected static pressure at the model is 

Pc * i3.^) 

where is obtained from the tunnel calibration, based upon contraction 
entry piezometer pressures. Figure 3.^ shows typical values of blockage 
ratio as a function of nominal jet velocity ratio R 

NOM * 

To complete the free stream data, tunnel density is obtained using p 
the measured tunnel temperature and the equation of state. Tunnel pressure 
ratio (to p^) , true speed and and blockage ratio (i.e. Ur / U„) are also 
calculated. ° 


3.3 Jet Hardware and Its Calibration 


Figure 3.5 shows a longitudinal section through the jet plenum and 
flow conditioning system with a 3-inch transition piece to a 1-inch Inclined 
jet fitted. For vert ica 1 jets , the transition piece is omitted. In the 
latter cases pi enum-to- jet pipe contraction ratios are 12.8 and 115 for the 
3-inch and 1-inch pipes respectively. Pre-straightening, from the supply 
cones at the base of the plenum (Figure 3.6) Is accomplished by a 2-inch 
thick honeycomb of i-inch cell size. When installing the inclined, 1-inch 
pipes the plenum position was adjusted to keep the jet exit location constant. 

Calibration runs, using a total pressure rake at the jet exits, showed 
a need to correct the 3-inch jet profile to match the more fully developed 
1-inch jet profile. Figure 3.7 shows an adjustable, profile-modification 
device used to increase the velocity deficit around the 3-inch pipe. A 
g^d (^atch was obtained, after two adjustments, for the configuration shown. 
The profiles themselves will be discussed In Section A. 


Jet calibration 

Pressure ratio, to corrected tunnel static, p^, was the primary jet 
control variable. Mass flow was also measured but this included roof 
boundary layer control air, during forward speed runs with the larger jets 
at high velocity ratio. The main calibrations were therefore made on the 
basis of jet p lenum-to-ex i t pressure ratio. 


A rake of total pressure probes at 0. 1-Inch spacing was employed to 
determine the thrust of the partially developed pipe flow profiles at jet 
exit. Data were obtained across two diameters at right angles and integra- 
tions were performed to determine mass flow, for comparison with direct 
measurements, and thrust. Checks on the thrust gave good agreement, for 
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partially developed pipe flow skin friction, with data quoted by Ower and 
Pankhurst 


The following calculation procedure was employed during data 
reduct ions : 



(3. S'! 


jet 


/(I 



M. ■') 
jet 


(3.6) 


{3 T. 
jet 


th 


M. . >R T. 
jet jet 


The mean jet velocity Vj was obtained durinq calibration from the 
measured mass flow in via: 


•J Aj 


A shape parameter, \^was defined as 


th 


(3.7) 


\ J . 


(3.^) 


L3. 10) 


\ is obtained during jet calibration and is used routinely to obtain the 
flow coefficient Cq using 


r i ^ i 

« Aj UcS Aj 

• J S 

■ .’c S ■ Aj 


(3. Ml 


where \ is a function of jet pressure ratio (see Figure 3.(3) 
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As net jet thrust could not be measured directly, and I 2 , defined via 
equations 3.12 and 3.13 below, were evaluated from calibration rake data: 


* 

m ■ 2”pj 


V3_ r dr - m 2-ipj 


^ r . r 
R 


^ dr- -d,2vjVjI, 


\l ^ ''m ''m 


(3.12) 


T ■ H2 Z’^pj 


J /Vj ■ -2 

— -I -r- d 7 T- “ n.,2:rp|V,I_ (3-13) 


''Jra ^ 


- R, *2--J-J*2 


T) and n, are calibration factors which allow for the differences, due to 
traverse^coverage and other reasons, between rake integrations and true 
values. Since m is known, we may find nj via (3.12), from 


m 



2^PjVjI, 


( 3 .IM 


The final step contains the main assumption of the present calibration, 
that n, and n, are equal; i.e. the calibration factors for mass flow and 
momentum flux^are the same. This assumption is supported by the fact that 
good pipe skin friction estimates were obtained from the present cali*- 
bration data. 


we may combine (3.11») and (3.15) and substitute for in (3.13) to 
obtain 


(3.15) 


T - 






mV 


mV I 


. ^2 

From (3.10) T - mAV^^ — 


I 


( 3 . 16 ) 


Finally, we may obtain the momentum flux coefficient Cp S/Aj from 
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mV. 




S 

A 


J 


N 


mV 


th 


7'c‘^c‘^J 


(3.17) 


>2 . «£ 


and (i.e. are shown for th'' three Individual configu- 


th 


rations in Figure 3.8, as a function of jet pressure ratio. It will be 
noticed that \ and Vg/V^h are much lower for the modified 3*inch jet than 
for the clean pipes. This is because of the extra pressure drop caused 
by the flow shaper in the former case. 


3.^ Jet Mass and Momentum Flux Coefficients at Forward Speed 

Figure 3.9 shows the nominal true velocity ratio, Rnqm* a 
function of jet pressure ratio, for typical runs at forward speed. 

Decreased pressure ratios, for 3-inch jet cases at « 8, reflect the 

fact that these tests were run at a lower mainstream speed than the others. 

In the case of the modified 3-inch and the 1-inch jets, R^qm was 
found to be a very close approx i mat i on to the effective true velocity ratio 
Vf/Uc- This was also true for the clean 3-inch jet, up to RnOM*^* 
thereafter Vf/Uj. increased faster than Rnqm. to 8.25 at Rnqm " 8.0. 

Differences between jets were noticeable for both the flow coefficient 
Cn (Figure 3.10) and the momentum flux coefficient C,; (Figure 3.H) when 
plotted as functions of the true velocity ratio Rnom* However, these are 
the true test values applicable to RnOM values used in comparisons of jet 
trajectory and vortex strength information. 

It has been found that the spread between large- and small -jet curves 
in Figures 3.10 and 3.11 is due predominantly to the fact that the large 
jets ran 10"'’-15°F cooler. The effects on density ratio are shown in Figure 
12. If RnOM 3 nd weighted by the density ratio, as in 

Figures 3.13 and 3.1^, a much improved correlation is obtained between the 
modified 3-Inch and the 1-inch jets. The differences which remain reflect 
firstly the fact that there are some slight Inconsistencies amongst the 
interim jet calibrations used for the various pipes during test runs. 
Secondly, C.^j, In particular, should not correlate on an RNOM^s^^^Fed basis 
for differing profile shapes: the distinctions between average, RMS. and 
jet transport velocities become Involved. 
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It may have been noticed that the calibrations and analyses described 
above were made with especial care. This was necessary to Insure that com- 
parisons between large and small jets will truly reflect tunnel effects 
rather than differences between the jets themselves. 


3.5 Wall Pressure and Jet Traverse Measurements 

Five sets of pressure orifices were installed, as indicated in 
Figure 3.15, in rows which extended for most of the test section length. 
Symmetry checks were possible using rows 3 and 5. However it was not 
possible at the time of the test to install orifices on the right hand 
wall, because it was all-glass. 

Figure 3.1 shows the rake of seven pitch-yaw, five-holed pressure 
probes used for jet flow traverses. This rake was mounted with the probes 
horizontal and the rake axis vertical for all tests. This limited in-jet 
traverses to locations no less than 6-diameters aft of the jet axis for 
the velocity ratios of interest. Other traverses, just aft of the jet 
figure 3.2), were restricted to a small region centered in the plane of the 
jet exit. 

The traverse measurements fell into two distinct sets. The in-plur 
measurements at X - 6D were for the determination of vortex strength and 
location in the plume. The remaining measurements, in the jet exit plane at 
X - 0 and 6D were for use in matching jet velocities during plume modeling. 
The two X-loca:ions may be regarded as representing "wing'= and "tail'' 
locat ions . 

Data were recorded and reduced using standard techniques, though 
second-stage analyses-invol ving stream function and circulation calculations 
were fairly elaborate. These will be described more fully in Section . 
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4.0 TEST RESULTS 


4.1 Similitude Between 1- and 3”inch Jets 


Jet pro filet? 

Figure 4.1 shows typical velocity profiles, measured during static cali- 
bration, for the unmodified and the modified 3-inch jets. For the 3-inch, 
unmodified configuration, the core flow profile was tilted. The axis of the 
tilt could be changed by biasing the supply (Figure 3.6) but it could not be 
eliminated. Nonetheless, the velocity in the central core was uniform to 
within about ±2% and the difference between the clean and matched profiles 
was sufficient to be useful for investigations into the effects of profile 
shape. 

Figure 4.2 gives a comparison between the modified 3“incl' jet profile 
and the 1" pipe profile to which it was matched. It is evidert that a good 
match was obtained. 

Botto pvcse:a't'o on Jot oylinJcvs 

The cylinder diameters and forward speeds used were suc'i that the Reynolds 
numbers for both pipes lay in the high-Cp, pre-transi t ional range. However, 
the possibility was recognized that finite length effects might reduce the 
drag coefficient of the 3-inch pipe (see Prandtl and Tietjens - Applied Hydro 
and Aeromechanics, p 97). So checks were made on cylinder b-jse pressures. 

Figure 4.3 shows the base pressure distributions down each pipe for 
three typical velocity ratios. Good matches are achieved between the 1-inch 
jet and the with-BLC, 3-inch jet data. The base pressures attain the two 
dimensional value about 5- to 7-diameters below the jet exit. The reason for 
the shift in the no-BLC, 3-inch jet cases is not understood. However, the 
magnitude of the shift is not enough for it to be of great concern. 

It appears that end-effects, on the 3-inch cylinder, do not cause serious 
flow changes, relative to the 1-inch data. Cyl inder- induced effects upon the 
jet development should therefore be the same in both cases. 

4.2 Traverse Results: Flow Distributions 

Heasurements were made on a sample basis in the longitudinal and hori- 
zontal planes and on a comprehensive basis in a transverse plane at x«60. 

The longitudinal and vertical runs were to explore the general nature of the 
flow. The transverse traverses provided vortex location and strength inform- 
ation for flow modeling purposes: the traverse boundaries were selected 
accordingly. Oownwash data, in the jet exit plane, intended for use in modeling, 

will be summarized below. 

Lo>uj : t :ui: I p lone (Y = 0^ 

An initial series of traverses, not mentioned previously, were made in 
the longitudinal center plane of the jet using a laser velocimeter. Anomalies 
found in these data, which were traced subsequently to a loose mirror in the 
LV system, made retesting necessary using the pressure probe system described 
above. The change restricted quite severely the choice of practical traverses. 
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Though the LV data were useless as 'hard' data, for finding cross-flow 
vorticity distributions, etc., they did give general guidance concerning jet- 
trajectory in side view. 

Figure ^.4 shows velocity vectors for the )-inch jet at RnOM“ 3. A good 
impression of the flow structure is obtained, both in the jet bending region 
and behind the cylinder. The flow vectors aft of the cylinder suggest clearly 
that the increase in cylinder base suction towards the jet exit (Figure A. 3), 
is associated with jet-induced vertical flow. It is noted that Figure 
shows only the upper four diameters of the 1-inch cylinder. 

Total velocity contour plots corresponding to Figure k.k were compared 
with published total pressure data at the same velocity ratio for a jet 
emergent from a flat plate. The comparisons made, at several R^OM values for 
the 1-inch jet, indicated that the trajectory shapes were similar and suggested 
that the form of the trajectory equations developed for jets-f rom-surfaces 
should be applicable to the present data. 

Horizontal plane (Z-0) 

Axial velocity distributions near the sides of the cylinder (Figure 
4.5(a)) show the expected, doublet-like character. The wake region fills 
and spreads quite rapidly opposite to the jet, which induces a strong down- 
wash component there (Figure 4.5(b)). The dov/nwash reaches a maximum about 
three diameters aft of the jet. The zero downwash contour at X=» 6-inches lies 
just to the outside of the plan-view of the vortex trajectory (see below). 

Figure 4.6(a) shows downwash at four locations in the jet exit plane 
as a function of RnOM* upper plots are for a location one diameter aft 

of the jet, representing the near field; the lower plots are at 6 diameters 
aft of the jet. The span locations were chosen so that comparable data were 
available for both 1- and 3-inch jets. Downwash values quoted are means of 
left- and right-side data in all cases. 

On a simplified basis of classical ang 1 e-of-attack correction theory, 
about one degree of downwash reduction should be experienced by the 3“it>ch 
jet at Rmom* 3- higher jet velocities three dimensional and impingement 
effects make such estimates inappropriate. Though Figure 4.6(a) exhibits 
significant data scatter, it is apparent that the anticipated downwash re- 
duction for the 3“inch jet is absent. In three of the four plots the 3"inch 
jet experiences more downwash, not less, than the 1-inch jet. 

Figure 4.6(b) shows similar data for the inclined and vertical 1-inch 
jets. The differences between jets are more marked, probably because of 
changes in vortex spacing (see below). 

Transverse plane (X»6D) 

The main 'production' measurements were made in transverse planes six 
diameters aft of the center of the jet exit for five configurations (1" 
vertical, two 3” vertical and two 1" inclined jets) at five velocity ratios 
(R«2,3,4,6 and 8). The primary flow measurements were total pressure , static 
pressure, axial velocity, vertical velocity and lateral velocity. Axial 
vorticity distributions (the curl of the transverse velocity field) and source 
strength distributions (the divergence) were derived. Streamlines were 


calculated from the axial vorticity distributions: this was a more con* 
venient procedure and gave more intelligible results than working directly 
with measured cross-flow velocit/ data. 

First-stage 'production' plots for the twenty-five configuration/velocity 
ratio combinations comprised: total pressure, vertical velocity and lateral 
vtlocity contours, cross-flow velocity vectors, vorticity and source strength 
contours . 

Only the more useful plots, of total pressure, vertical velocity, lateral 
velocity and streamlines are presented in this report. These are consolidated 
to show the effects of jet inclination in Appendix A and to show the effects 
of profile shape and jet si?e in Appendix B. Summary data, taken from these 
plots, will be presented in subsection ^.3. 

Figure ^.7 repeats selected data from Appendix B and shows the effect 
of jet si2e upon cross-flow streamline patterns as Rnqm increased. The 
effects of profile shape were found to be secondary. Length scales in 
Figure 1*.7 are measured in diameters. The tunnel roof location, indicated 
by the 'hatched* regions, consequently appears closer to the jet exit for 
the 3-inch jet cases. Negative stream function contours arc denoted by broken 
lines and the contour values are listed to the left of each plot, starting 
with the innermost negative contour. The integers to the upper left of each 
plot are run numbers. (Test MTF68) . 

At R ■ 2 and 3 (Figure 4.7(a) and (b)) the vortex penetration into the 
flow has apparently not been affected by the proximity of the tunnel roof 
to the plume for the 3-inch jet though there is somewhat less vertical elong- 
ation of the vortex streamlines. At R * 4 (Figure 4.7(c)), both a reduction 
in vortex penetration and a marked flattening of the st'eaml ines is evident 
for the 3-inch jet. Another indication of roof-induced interference is the 
reduced magnitudes of the stream function values at the vortex centers - a 
consecuence of a reduction in vertical velocity between them. 

In the R«6 and 8 cases (Figures 4.7(d) and (e)) the 3" jet impinges 
strongly upon the tunnel roof and the flattening and the other effects noted 
above are very pronounced At R*8 there was significant flow unsteadiness 
not only for the 3-inch but also for the 1-inch jet. Application of tangential 
blowing of *’he roof reduced the unsteadiness, though a tendency for the 3-inch 
jet to bend sideways may still be seen. 

The data of Figure 4.7, taken with data in Appendixes A and B and else- 
where, indicate that the structures of the 1-inch and modified 3-inch jets 
are qualitatively similar at low values of jet -to-ma i nst ream velocity ratio. 

It may also be inferred from experience with the 3-inch pipe, that streamlines 
for the !-inch pipe at R«9 are sufficiently far from the tunnel roof, as in 
Figure 4.7(e), so that tunne I - i nduced local distortion may be assumed to be 
neg 1 ig i b 1 e . 

4.3 Traverse Results: Data Summary 

-it X ■ 6D 

Having calculated the streamline patterns, vortex centers, defined as 
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maximum- or minimum-; points, may be found quite '■eadily. Figure 4.8(a^ 
shows that vortex penetration into the flow, for the 1-inch jet, is very 
nearly proportional to RnOM- Comparison with Weston's data (Ref. shows 
less penetration for his data, which are for a jet emergen.* from a plane; 
this seems reasonable. However, it is pointed out that there are differences 
between the present analysis and Weston's in the definition of vortex center. 

The two sets of 3-inch jet data differ only slightly from each other, 
mainly in vortex span, confirming the insensitivity to profile shape, 
mentioned earlier. The jet penetration curve for the 3*inch jets diverges 
from the I-inch data at about Rnqm * 3 as the tunnel roof causes increased jet 
turning. It appears that the 3“inch cl rve will asymptote to about 0.7D below 
the roof at high jet velocity ratios. 

Figure 4.8(b) shows that forward i nc 1 i na t i on of the I-inch jet causes 
increased spreading of the vortex pair, but reduced penetration into the flow. 
It will be seen later that the increased spreading is accompanied by reduced 
strength, relative to the vertical jet, at high Rnjqm- f^NOM * . 3 and 4 

vortex penetration for aft-inclined and vertical jets is almost the sane. 
Thereafter the aft-inclined and forward - i nc 1 i ned jet data are tbe more similar. 
Vortex spreading is less for the aft-inclined jet than for the others. 

dt X ■ 6D 

Aftfr considering a number of alternatives, a method of determining 
vortex strength was selected based upon line integrals of velocity around the 
streaml ne patterns derived from crossflow vorticity data. This is consistent 
with the vortex center determination procedure, just described, and also per- 
mits an economical description of diffuse vortex cores. By plotting vortex 
strength as a function of the stream function (Figure 4.9) a family of similar 
curves is obtained with R^qm ® parameter . Intercepts on the horizontal 
axis are vortex center locations, in terns of ,. The inaximum circulation 
strength is usually associated with the zero or a nearby streamline and is of 
prime interest in the present work because tunnel effects relate to the far- 
field. Central streamlines sometimes intersect the traverse boundary, which 
is used to complete the integration circuit in such cases. This provides 
the opportunity to calculate both clockwise and anticlockwise circuits for 
the same central streamline and can give rise to the overlap region (both 
open and filled points at given v) which is evident in Figure 4.9 at Rf 40 M =* 6. 
The occurrence of such a discrepancy is an indication of net circulation 
around the outer boundary of the traverse. Beth ^rnax po'hts will be shown 
subsequent 1 y . 

Figure 4.10 shows maximum values of vortex strength obtained from Figure 
4.9 and other, similar plots. It should be noted that the data represent 
the axial component of the vortex strength, not the total. Despite differinu 
exit conditions and differing methods for finding vortex strength, it is 
found that the present results ar* quite similar to corresponding data, 
derived from Ref. 4, for jet -f rom-surface cases. Since vortex strength is 
greater in the present case, and spacing is less, their product - which is 
proportional to lift - tends to correlate better. Two apparently 'bad* 
points, at R ■ 8 in Figure 4 10(a), were derived from a . curve which 

had an Inverted peak, uni ike any ot*'er data. This probably reflects an ou . - 
of-range condition for the pitch prole, so these points have been ignored. 
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The 3”inch jet circulations in Figure A-lO(a) agree well with the 1-inch 
data at ^ somewhat lower at Rnqm “ 3 and k. As impingement 

intensifies, a circulation plateau is reached at about 2.5 U„D. 

Complementary behavior between vortex strength and span is again ob- 
served for inclined jet cases (Figures 10(b) and ^.8(b)). After behaving 
similarly to the vertical jet at R“2, 3 and k, vortices in the swept-forward 
jet display 'ncreased span but decreased strength. Conversely, vortices in 
the swept-afr jet plume are less widely spaced but stronger than for the 
vertical jet. 

Vortex-lift correlation 

The vertical component of jet reaction lift is given by 

where 6 Is jet inclination to the vertical 


and 


2C 


2l 


ipU 2A 


Poo [Uoo j 


cos6 * R^cosS 


(^. 1 ) 


Vor tex- 1 i f t , from the wake measurements, is given by 

1-2 “ ^®U»*max“'^ 

and 

''^2 = IpU^Aj ■ ° U„D D 

It is evident from Figure . 1 1 that, for the 1-inch vertical jet, 

C 2 . • C 22 to a reasonable approximation, i.e. all of the jet reaction lift 
i s "represented by the vortex pair. For the inclined jets C 22 > Cj,. The 
Implication ray be that the vortex pair represents total jet thrust rather 
than the lift component. Removal of the cos5 factor from C 2 ^ improves the 
correlation in Figure ^.11 but does not collapse the data entirely. 


The above result is significant in relation to jet- I n-crossf low modeling 
because it provides a rationale by which vortex strength may be estimated if 
the spacing is known. 



Tunnel Surface Pressures 


ru*u' .:r .h,. 


S- ^ 


Figure h 12(a) shows tunnel surface pressures with the 3-*nch jet pipe 
installed but with both the blowing air and the ground BLC turned oft. The 
upper part of the figure is an empty-tunnel -plus pipe case and shows suet. on 
peaks, corresponding to the maximum cy I i nder/separat i on bubble diameter, 
followed by a pressure rise towards an asymptote proportional to wake dis- 
placement thickness. On installing the traverse gear (lower plot) this 
pressure recovery is largely lost and further flow acceleration occurs 
opoo^Ue t^he traverse' geir location. Figure V12 M shows the prev.ous 

distribution is nedified when Aground* BLC is app..ed. Figure 14.12(c) shows 
corresponding data for the vertical 1-inch Jet pipe. Here, the effect ot 
the pipe itself is quite small, as expected, and the traverse gear signature 
predominates. The effects of sweeping the jet pipe 30 -degrees forward or 
aft are found to be smal I . 

In subsequent figures the appropriate jet-off datum values, depicted in 
Figure k 12 have been removed in most cases. Certain exceptions occur m 
strong impino-ment regions where greater- than-ma i ns t ream total pressures 
invalidate the »^(T^Cp) superpos i t i on technique which was used 
are recognized readily since Cp : 1 at impingement. 


These rcoions 






This jet intpioqes at about RhOH " '’“h ‘ '*• 

ment is present (Figure It. 12) and there is evidence of a separation vor u k 

ahead of the X - 0 station in the three no-BLC cases shown. Rpnlvmo BLC 
(trianoles) destroys this vortex and may reduce the impingement pressure. 
Comparisons between sets of octagonal points show good repeatab i 1 . ty between 
runs and comparison of these points with the 'plus' points suggests that the 
effects of jet profile shape are not very great. 

Figure li.l'* shows similar. no-BLC data for various RnOU VALUED- 

R„pm - C and 8 the vortex suction peaks rove forward and becomes " 

i 5 °Slso evident that impingement Cp's become very high. There was cons ii.erable 
unsteadiness in the flow at Rnom ‘ caused scatter m the data. 


3- : 


tu*i>u'2 .vvi.’r prrssurre 


At Runn • 2 Figure A. 15(a) shows that the wall pressures are blockaoe- 
dominated: °there is a continuous acceleration along the 

imposed upon this, comparison between rails 1 and 3 or 5 hut then 

anticipated difference signature associated with positive jet l.f but t . 
downstream of about X - 0.2. a reversal which corresponds to negative jet 
lift Thii may be due to the fact that jet blockage effects are offset 
towards the ’'ground” (i.e. the tunnel roof). Comparison between the upper 
and lower plots in Figure 4.15(a) shows that ground BLC has little effect ut 

R - 2. 
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^NQH * ^ (Figure 15(b)), spreading between rails 1. 2 and 3 becomes 
more noticeable as lift effects increase. Application of 'ground' 8LC (lower 
plot) reduces the spread significantly. This may be due to entrainment into 
the BLC jet slieet, which will have an ac~eleratinq effect upon the flow at 
rails 3 and 5. 

B'' ^NOM “ ^ (Figure 4.15(c)) tunnel flow breakdown well advanced 
and the locus of the characteristic impingement "footprint" and the vortex 
ahead of it may be deduced quite readily. Since the vortex ts undoubtedly 
skewed quite strongly, the X/B-length scales in the pressure plots may exag- 
gerate the vortex size. Rail 3 and 5 data, at Rnqm " ^ suggest that sepa- 
ration, vortex center and reattachment are at X/B ® -0,20, + 0.15 and * 0.70 
approximately . At mid-sidewall (Rail 2 - triangles) the same features may 
be identified, shifted downstream, but only the vortex peak location can be 
estimated, at about X/B » 1.0 (upper plot. Figure 4.15(c)). This was 
changed somewhat by 'ground' BLC (lower plot) but otherwise the BLC produced 
no great effect. Some differences are observed between right and left sides 
(Rails 3 and 5) at the two highest R^qm values (Figures 4.15(c) and (d)). 

This is most likely caused by the jet bending sideways (see also Fiqures 
4.7(d) and U)). 

of pvcf^'Zr shive 

Comparisons between corresponding plots in Figure 4,15. which is for 
'square' profiles, and Figure 4.16 ('pipe' profiles) show that differences 
due to profile shape are almost negligible. 


Figures 4.17(a) and (b) comprise four pressure plots for Rnqm values 
of 3. 4 and 6, respectively. Curves for several jet inclination ang’es are 
shown on each. The overall magnitude of the signatures are small prior to 
impingement but the general trends are the same as those described above. 

As might be suspected, the f orward- 1 nc 1 I ned jet (5 ■ -30“deqrees) produces 
greater effects at the walls than the vertical jet and the af <■- i nc 1 i ned jet 
produces smaller effects. 

Figures 4.l8 and 4,19 include data for other rails at R^PM ^ 
respectively. In the “ 4 case the magnitudes of the signature are 

starting to become comparable with measurement accuracy and flow unsteadiness 
effects. Nonetheless, the anticipated trends are present. The Rnqm “ ® 
cases (Figure 4.19) f’^ay involve impingement and both lift and blockage effects 
are starting to become noticeable. 
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5.0 THE 'VSD‘ now MODEL AND IT'S DEVELOPMENT 

5,1 Int roduct i on 

The present aim is to identify or develop a realistic theoretical 
model of a jet* in-crossflow whi ch is suitable for tunnel interference 
estimation. The model will be applied with the standard wall pressure signa- 
ture procedures, firstly to extract jet effects from powered-mode 1 -s i gnatures 
and then to provide jet-related interference velocities for addition later 
to those determined for other parts of the test model. 

Any flow model used ^or these purposes must be reasonably simpler this 
eliminates the finite difference and the detailed vortex lattice techniques 
used for many near-field studies It was thought at first that an extension 
of Heyson's work, probably towards the Fearn, curved trajectory model would 
be adequate. However, both were found to have serious shortcomings and 
considerable development work was needed to produce a theoretical model 
which was reasonably simole, yet reptoduced observed flow features oroperlv. 
The new flow model, desiqrated the 'VSD' model (Vor tex , Source , Doubletl will 
be described in subsection 5.2. 


The terms of reference for the present work exclude impingement cases. 
However, cases in which the theoretical plume strikes a tunnel surface 
cannot just be ignored. Appropriate geometric changes have therefore been 
made, for impinging cases, which parallel Hevson's treatment. Though this 
approach appears reasonable, detailed studies have not been made of the 
impingement region itself or of ways to model it properly. Results obtained 
here for impinging cases must therefore be considered speculative. 


5.2 Description of the 'VSD' Model 


The data found most useful for modeling purposes were measured vortex 
strengths and locations - augmented by the 'Fearn' vortex model - and tunnel 
wall pressures for the 3*Irtch jet. 

As the work progressed, it became apparent that twin vortex models 
such as Fearn's contributed only weakly to the wall pressure signatures: 
Sources and/or curved lines of axially-directed doublets were needed to 
match the observed tunnel blockage effects. As the main body of the jet 
penetrates further than the vortices, the source/doublet lines were given 
the greater penetration into the flow^as sketched in Figure 5.1 (a) and (b). 
On the basis of measured data at X/D “ 6 and 2 < R < 8, the 'Fearn' and the 
'Williams and Wood' trajectories respectively were selected for vortex and 
source/doublet elements (see Figure 5.1(c)). 

The line vortex and line source/doublet trajectories are defined for 
the 'VSD' model via the equations; 


/■ 
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7/D - 0.35 


(5.1) 


for the vortex pair 


Z/D - 0.758 (X/D) ^1.000 2) 

for the source/doublet pair. 

and 

Y/D-O.O/b^^ (X/D) f^l.OOO 

for both vortex and source/doublet 
pa i r 5 , 

I’v'rrt'j: strrK.jth 

Since the Fearn equation for vortex strength agreed well with the 
present data, this Is used directly in the 'VSD' model, i.e. 


r 

U.vD 


0.600 


R-' 

(xTd) 


(1 


-0.035(X/D)-\ 
e ) 


(5.M 


The source system is selected so as to model the physics of the 
observed wakes. The system described is the original one, derived as stated 
no adjustments to constants were used to improve the match to the wall 
pressure signature. 

For elements of the initial part o' the plume, the data show that an 
almost para I le I -s i ded wake is required (e.q. Figure A1, all R) . The total 
volume Qi 2 emanating from a line source between points 1 and 2 in the 
initial plume may be written 

Q|2 • .^Y (Z 2 - Zi) (5.5) 

where AY is the far-wake displacement width. For the initial plume, a 
value AY • D (i.e. ■ 1) appears to be appropriate. For a plume 

developing without a crossflow, entrainment adds to the mass flow in the 
jet at about a 23v rate. An improved approx imat ion for AY is therefore 

AY - D .'(1+0.23 Z/d) 

Recognizing that mixing depends upon path length, S, the flow physics 
Is better represented by 

AY - 0 r^o.23 S/d) 

Substituting into (5.5) we obtain 



The local source strength, per unit length along the plun,e. is then given 
bv 

Ut^D (22 ’ 


Qi2 


12 


.^S 


12 




r^O.23 S 12 /D 


(5.6) 


Doublet strength 

In addition to the sources, a doublet system was introduced in 

response to lol irbCka^^^^^^ present, which 

case, local peaks made it upstream-directed axial 

requires either a ji^^^ter at exit and grow at 

?hfl:m:-rale“t f7r system. The-axial doublet strength distribu- 

tion is then given by t 

(5.7) 


12 


-j U.x,D‘ I' (1+0.23 Si 2 ^d} 


No vertical doublets are required because lift 

system. In evaluating (5.6) and (5.7) the approx 1 mat . on 5,2 X,j 


5.3 Tunnel Surface Pressures 


tne 


nooe . 


D’/^itbes 1 s' 

The development of the VSD model and the i' 

are best ' O^'e^^no’’ imoi ngement ) will be used to 

The case with the 3 "'nch jet at R - ^ U-e. no y 

demonstrate the matching procedure. 

r \ u c that Fearn's vortex model causes almost negligible 
Figure 5.3(a) shows tha 3 „<j represents rail 6 con- 

effects at the tunnel s • dewa J the vortex. drag 

ditlons poorly. It may be nri<;<iures very little. A separate 

implied by the Fearn model affect the doublets in a flow 

e-iimate was made of vortex drag and used to s e „ 3,3 

model like Heysons (see Appendix C). The results IFigure 5. at 
similar to those for the Fearn model. 

in view of this failure of ‘ we^e“tl!est igated 

pressures, the properties of ^ ^ conslant-strength 

when laid out along a curved jet ^ Houblets sized to give far 

sources or J (Equations (5.6 and (5.7)) 

wake width D and cyl.nder ^ ° 5 I ?he plus symbols, for 

^^rder:orrr;Fq::^:^s: 6 :t::’;:;"^j^^^^^ 

?:rL::;'rr nril^terr-unsrele-rra^-riimnar to those for the 
doub lets. 


19 



f 

\ 

\ 


Comparison of Figure 5.^ results with experimental data (Figure-^ 

5.2 and 5.3) showed clearly that at least sources are needed to complete the 
flow model. On the basis of wall pressure comparisons, there is a temptation 
at this point to omit the vortex system. However, this would be incorrect 
because the associated upwash inteference would be lost. 

Figure 5.5 shows measured wall pressures and predictions from a vortex 
pair/split line source model, with constant strength sources. The comparisons 
are encouraging. However, Figure 5-6 shows that some improvement is possible 
bv grading the source strengths as indicated earlier. Because of their sim- 
ilar properties, it is evidently possible either to Increase the source 
strengths or introduce upstream-directed doublets to improve the correlation. 
The latter choice was made on physical grounds. 

Figure 5.7 shows that the introduction of doublets, which completes 
t^'e VSO model, generally improves the wall pressure correlations. Further 
fine tuning is obviously possible but was not considered worthwhile on the 
basis of one flow case. Comparison of Figure 5.8 with Figure 5-7 shows 
that the changes in wall pressures on removing the vortex system from the VSD 
model are almost negligible. However, the full VSO model was retained for 
the reasons indicated earlier. 


Figures 5.9(a) and (b) show wall pressures for the 1-inch jet at R = 2 
and R = 4, respectively. As the measured pressures are of very low level 
and are a residue which remains after removing the (larger) traverse gear 
effects, good comparisons are unlikely. With a single exception, however, 
the predictions appear qua 1 i tat i ve 1 y correct. The exception concerns Rail 
4 at R ■ 4 (Figure 5.9(b), which has obvious problems at lar e X/B. Figures 
5.10(a) and (b) show VSD and SO results respectively for an impingement 
case: the 3“inch jet at R ■ 4. Rails 1 and 2 correlate quite well but the 
Rail 3 ^nd 4 results show that the impingement flow is not well represented 
by the V r 3 mode 1 . 

flex: effeets 

In converting flo\v velocities generated by the theoretical models to 
pressures, it is necessary to use Bernoulli's equation and assume that the 
flow next to the wall has mainstream total pressure. This not so in the 
jet impingement region, as evidenced by experimental points near X/B * 0 ^or 
Rail 4 In Figure 5-10. It is apparent from the experimental curves that a 
stagnation pressure of about three times mainstream is present here. It 
appears highly likelv that the Rail 3 pressures are also directly affected 
by impingement total pressure effects, which are not modeled by the present 
scheme . 


5.4 Interference Velocities 

The tunnel interference flow, at the model centerline, comprises 
velocities due to the image system of the model and it's wake as it exists 
In the tunnel and velocities associated with any redirection of the model 
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wake caused bv tunnel surfaces. The second effect is usually '^nored. 

ptrurH.II inustrates this -^akdown for an i.o n,,n, ,e^ 

irrShfehut‘?h:";™:d 're^:n!'1ho:^h"hesrefrect^l!rL'e^ 

model the impingement region itself. 

The remainder of this subsection will comprise a description, within 
the above Tramework, of interference predictions first for the 't;- 
then for the 3- inch jet using the VSD flow model. Compar.sons w,ll then 
made with similar predictions using other methods. 

(Figio'C 5.12) 

The R - 2 li and. to a lame degree the R » 6 cases for the 1-inch 
^ in the present work because impingement is not 

e c mjll a. connected 

ment and redirection (see above) do not aris . y 

iiitiiv. ai o 0 = 7 and k Increase with X/B in mt 

:trt?r:ar:rasv:™:4:r:priately far d,«nstream (figures 3.l2(al. 

(b) . and (d) ) . 

At R - 6. non-planar effects become important for j*'; 

s:;:r;w;s^;i:r;rrar:r*'i:c:: ::;^i:j!auark'i™ts ingures s.^cai 

and (b) ) . 

Figure 5.12(c) shows the '"crements associated ' 

to be added to the previous, image- .nduced ' ^ impingement 

involves swinging the plume throug ,3 3 , ,.05B and the 

--."X^irii-de^r!^!; The effect of mdi rec. ion^i s m increase the 

lratlinr:"rex^:irr^:;res^?::"rrl^ 

redi rect ion. 

The 'bottom line' for angle-of-attack is shown In f' 3 uce 5.12(d). At 

;rn;rp;:::ri:a;: rff:-::"! i:cr:;rerfa:::r;h:rR-sguared,the 

rate for planar assumptions. 

Jft (Figure' 5.1^' 

Fnr this let only the R » 2 case is non- i mp i nq i og : most of the 
° ‘ 1 : n case for the renaming, impinging case . 

previous comments apply m thi. ca . further discussed 

the previous cautionary remarks apply. The ° 

because of some interesting findings concerning the redirection term. 
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In the absence of the along- surface elements (Figure 5.13(a)) the 
interference effects are quite insensitive to R. Once impingement occurs 
(R ■ 4, 6 and 8), the VSD model comprises a curved cylinder joining the Jet 
exit to the Impingement point, which has limited movement over this R-range, 
Blockage effects are therefore almost the same at R - and B (Figure 5.13(a)) 

The movement of the angle of attack peak follows that of the impingement point 
as it progresses forward with increasing R. 

The introduction of a long-surface elements (Figure 5.13(b)) and their 
subsequent redirection (Figures 5.13(c) and (d)) have major, but opposing, 
effects upon the corrections. The two corrections illustrate the importance 
of redirection vividly. Figure 5.13(h) shows apparent interference calculated 
for a ' ki nked-plume ' case (Figure 5.H, sketch (I)). Here the plume is held 
next to the tunnel surface artificially and high ang 1 e-of-at tack increments 
result (compare Figures 5.13(3) 3 nd (b)). Since the true plume position is 
much further from the tunnel centerline, most of the spurious t ra i 1 er- i nduced 
effect in Figure 5.13(b) is removed when the plume redirection term (Figure 
5.13(c)) is applied. As a result, the magnitude of the net interference 
(Figure 5.13(d)) is several times less than some of it's constituents. This 
is not a reflection of any real effect but rather a demonstration that the 
intermediate, along surface model of Figure 5.13(b) (or Figure 5.11. sketch 
(it)) is inappropriate. .Methods which fail to pick up the redirection term 
must be considered suspect or, at best, incomplete. 


5.5 Comparisons Between Methods 


Comparisons will be made firstly on the basis of effective source and 
doublet distributions (Figure 5.1^), then in terms of interference predicted 
by the VSD, Heyson and other methods (Figures 5.15 and 5.1b) and finally In 
relation to results calculated directly from wall pressures using the 
program from Part I of the present report (Figure 5.17). Further comments 
on Heyson's method, as interpreted herein, are given in Appendix C. 


t' jruz distiKDUtioy'.s 


To illustrate the nature of the 'VSD' and Heyson flow models. Figure 
5.1k(a) and (b) expresses doublet effects in terms of equivalent circular 
cylinder diameter, plotted with the trajectory in the upper parts of the 
figures and expresses source effects In terms of displacement. 

For the 'VSD' model (Figure 5.1^(a)),the effective cylinder diameter 
(upper plot) is essentially that of the jet: the spreading term is weak. In 
the lower plot, the full line shows continuously increasing mass flow due 
to the (explicit) line sources. The broken line shows the implied mass flow 
due to the point source effects at the joints between doublet lli.es (see 
also Appendix C). For the VSD model, this effect is related only to drag. 


The total source effect is shown by the chained line. It is evident 
that a rapid increase occurs in the first few diameters to a level somewhat 
exceeding the jet volume flow. Thereafter the rate of increase declines as 
the jet bends. 
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In Figure 5.1^(b) a comparable Heyson case is Illustrated for the same 
lift and for the same effective cylinder diameter. The source effects are 
much weaker because only doublet“related sources are present. The drag 
related source effect. In the vicinity of the model is weakened by a sink- 
effect at the jet exit caused by the lifting doublet system. 

Despite having forced the Heyson model to match the VSD cylinder 
diameter. It is evident that it suffers from serious shortcomings regarding 
source effects and their distribution. 

5 r i n t evfc 

Figures 5.14 and 5.15 show blockage and angle-of-attack Interference 
increments respectively for the jet at various values of R. The 

' Fearn ' -vortex model and the 'SD' source-doublet model combine to form the 
'VSD' model. Comparison is also made with results from the present inter- 
pretation of Heyson's method. 

Figure 5.15(a) shows u-component blockage interference, due to the 
tunnel image sets, for several theoretical jet- in-crossf low models. At 
R « 2 and R » 4, the interference predicted by the present method Is an 
order of magnitude greater than that of Heyson. This reflects the fact that, 
be’ng largely drag-ba<^ed for blockage the Heyson results have strong R- 
dependence. The SD and VSD present results, on the other hand, depend 
heavily upon measured wall pressures which have much less R dependence. The 
R 2 and R * 4 results are not affected greatly by the redirection term so 
the results shown In Figures 5.15(a) and (b) are almost the same. 

In Figure 5.15(b) it is evident that the vortex 'cross' effects upon 
blockage become more significant as R Increases. With no red I rect Ion 
term (Figure 5.15(a)) this effect is very strong, but spurious. However, 
the Figure 5.15(a) result at hiah-R shows how Important it is to include 
the redirection term for vortex models. 


The K * 4 case in Figure 5.16(a) or (b) characterizes tne various 
methods quite well. Though the 5D model has no vertical doublets, the 
angle-of-attack curve is very similar in form and magnitude to that of 
Heyson, except that the SD Impingement point is located (more correctly) 
further forward. 

As a result of vortex inclination near the jet exit and in the early 
plume, the vortex contribution (I.e, the Fearn result) is shifted aft at 
medium R-values. This causes an angle-of-attack plateau to occur where 
negative dAa/dx for the SD model and positive dAa/dx for the vortex models 
are about equal. This local detail Is undoubtedly very sensitive to modeling 
assumptions and should not be taken too seriously. 

As R increases, the 'ortex contribution (Fearn) to anq 1 e-of-at tack 
dominates increasingly, pa ticularly prior to redirection (Figure 5.16(a)). 

It is interesting to note that the magnitudes of the image-effect corrections 
are comparable with corrections calculated on a simple. .'(S/OC^ basis. 

The fact that redirection reduces the ang 1 e-of-at tach correction bv a factor 
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of three or four (Figures 5.16(a) and (b) . lower parts) may ^e 

preference of some V/STOL experimentalists not to apply (, mage-based) 
tunnel corrections. 

Direct use of Doll pressures and the Fart I metnod 

It is obviously possible to ignore the non-planar aspects of the jet- 
in-croLflow problem and use the Part I program directly to obtain a rough 
estimate of Interference. If the jet is small enough or if R is low enough 
this should produce acceptable results. Figure 5.17(a) confirms this 
one^noh jet tests at low R values. The blockage leve ,s ^ 
at most) and wall pressure signals are subject to scatter. Nonethele. s 
'VSD' model predictions (crosses) are in quite good agreement «"h the 
resuits from the matrix method of Part I. This agreement also t^at 

the 'VSD* method operates properly as planar conditions are approac 

Fiqure 5.17(b) shows the corresponding results for the three-iich Jet. 
Here the blockage is an order of magnitude greater. These comparisons are 
particularly Interesting because, though the 

at R-2, by the non-planar effects have become significant. ’ 

•VSD’ treatment gives a blockage curve which levels out as the je^ be 
(crosses). However, for the early part of the characteristic ^ 

jet penetration is incomplete, the results from the planar and non-planar 
predictions are in remarkable agreement. 

It is tempting to conclude that, if impingement is absent, direct use 

of the Part I method will produce good of thts 

probably true for blockage at the model position, though not a-t 
However; serious probler.is can be anticipated in est imat i ng anc 1 e-of attac 
interference on the basis of the Part I method if used with an influence 
matrix for center-tunnel elements. Runs paralleling those for 5- 7 

confirmed that the latter procedure yields entirely 5 ° 

attack estimates. This occurs because (as 

effects are almost "invisible” at the tunnel wall. The Part I method 
would respond predominantly to the "cross” effects of offset blockage an 
return corresponding, spurious ang I e-of-at tack predictions. 
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6.0 CONCLUSIONS 


6.1 Scope of the Present Work 

The aim of the work described, to adapt or develop a jet- i n-crossf low 
model for estimating wind tunnel interference, has been met successfully. 

A study of measured and predicted wind tunnel surface pressure signatures 
shows that the two most likely existing candidates - the models of Fearn and 
of Heysor - fall seriously short of representing the flow physics properly 
for the present application. A new, vortex-source-doublet (VSD) model was 
therefore developed with significantly improved flow physics. This incorporates 

o a curved, varying-strength vortex pair based upon the 
'Fearn' near-field jet-in-crossflow model. 

o a pair of curved, varying-strength source-doublet lines 
extending beyond the vortex pair but with the same hori- 
zontal development. 

o source-doublet strengths based upon viscous wake measurements^ 
jet cylinder diameter, and jet growth considerations. 

The 'VSD' flow model is employed in a pre-processor program which 
removes jet effects from whole-model wal 1 pressure signatures and provides 
jet-induced tunnel interference data. The revised wall pressure signatures 
and the interference data become input to the subsequent "rest-of -model " 
interference calculation described in Part I of this report. 

The above development was backed experimentally (see 6.2, below) and 
by a number of theoretical studies (see 6.3, below). It is not easy to val- 
idate methods such as the present one: 6.4, below, deals with this topic. 
Subsection 6.5 considers the special topic of impingement cases. 


6.2 Experimental Studies 

Though a great volume of jet-f rom-surface data exists, its relevence 
to the tunnel Interference problem was in doubt: a jet-f rom-cyl inder repre- 
sents a typical V/STOL configuration better. Experiments were therefore 
carried out on jets emergent from chimney-like cylinders mounted on the 
tunnel floor. The jets, o^ 1- and 3-inch diameter, emerged from 15-inch 
high cylinders mounted in the Lockheed-Georgia 30- by 43-inch wind tunnel. 

Three component velocity measurements were made in the planes X * 60 
and Z » 0 over a range of jet velocity ratios, z s R s 8. Two jet profile 
shapes were investigated for the (vertical) 3~inch jet and three jet angles 
were investigated for the 1-inch case (30-degrees forward, vertical and 30 - 
degrees aft). Especial care was taken in calibrating the jets. The experi- 
ments lead to the following main conclusions: 

The jet flow measurements at X s 6D showed that jet-from- 
cylinder vortex trajectories and strengths, for non- i mp i nq i ng 
cases, were in general agreement with jet-f rom-surface data. 
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o Results for square; and pipe jet exit velocity profiles, were 
essentially indistinguishable at the same Cu. 

o Comparisons of 1* and 3“inch jet data showed that vortex 

penetration (in diameters) into the flow was almost unaltered, 
during tests at Increasing velocity ratio R, until 'hard' 
impingement occurred. 

o Careful integrations of the X «= 6D cross flow data showed that 

the product of trailing vortex axial strength and lateral spacing 
closely approximated the jet lift. This demonstration Is believed 
to be a 'first*. 

o The above three facts were helpful when developing the flow model. 

Tunnel wall pressure 'signatures' were measured along five longitudinal 
'rails’ for all of the above conditions. The 3'irtch jet data provided the 
primary checks during VSD-model development. 

6.3 Theoretical Studies 

7he pZjiKar models 

Studies with the completed 'VSD' flow model revealed the following; 

o The contribution to the wall pressure signatures of the vertex 
portion of the VSD model is almost negligible: the lifting 
system is essentially 'invisible' to the tunnel walls. 

o Ang Ic-of-attack corrections due to the jet- i n-crossf low must 

therefore be estimated entirely on the basis of the jet-in-cross- 
flow theoretical model. 

Demonstration runs were made by submitting jet- in-crossf low wall 
pressure signatures to the (planar) Part I analysis program. These revealed 
the fol lowing : 

o At low-R for the 1-Inch jet, the VSD and planar programs gave 
blockage estimates which agreed well with each other, as should 
be expected. 

o For the 3“irtch jet the agreement was good in the vicinitv of the 
jet but, aft of the jet, a continuously rising blockage correction 
for the planar analysis was replaced by asymptotic behavior for 
the full, VSD model. The latter Is more credible on the basis of 
previous experience. 

o Predicted ang 1 e-of-at tack corrections, using the olanar method, 
were entirely erroneous. This was because the method interpreted 
offset blockage interference in terms of angl e-of-at tack , 

These results of the theoretical studies prove that the pre-processor 
approach to handling jet- i n-crossf low Interference is not just expeditious- 
i t i s essen t i a 1 . -yc 



Both Heyson's proqran itself and a new, equivalent program which 
extends it for wall pressure calculations were implemented. At a qualitative 
level the studies were fairly successful. Source*like properties of Ime 
doublets in Heyson's model were identified and Interpreted and some inherent 
shortcominqs, relative to the present, flow-based m.odel . were uncovered. 
Quantitatively, the present studies with Heyson's model were unsat • s^v i ng . 
largely because no balance data were available. The problem was compounded 
by difficulties in deciding upon ground rules for comparing force-based ana 
flow-based methods. 


6 . 1* Va ’ I da t i on 

A direct validation of a jet- i n-c rossf low interference moc'e 1 is .ikely^ 
to be difficult because of the presence of other components for any practice, 
jet-powered configuration. The V$0 model. In contrast to others reviewed, 
is Supported by the following facts: 

o V^all pressures are matched quite well for the 3“ inch jet. 

O For the 1-inch jet, wall pressure matching Is qualitatively 

correct and of the right magnitude: this Is the most that can 
be expected of the data. 

o It follows that far-field predictions by the method.and h^nce 
tunnel blockane nredicticns, are of tne correct order, 

o Oownwash differences, between Z = 0 data fc ' 1- and 3“if^ch jets, 

were too small to be of practical use in validating anq . e-o . -at tack 
inteference estimates, 

o Angle-of-attack Interference estimates, therefore , res t heavIW 
upon ream's vortex model, which is Incorporated in the 'VSD' 
nodel The Fearn model is supported by extensive test data 
gathered by the original authors and by more limited experiments 
described herein, 

o Comparisons with the planar method (see above), show that the 
VSD method performs properly as this limit is approached. 


6.5 Cases with Jet Impingement 
The experimenta* studies included a significant 


expc r I L a ' 

cases. Here, floor tangential blowinq was applied as 
the vortex induced suct'on peak ahead of impingement, 
a'-knowledqed Impingement but did not maUe any attempt - 

except with regard to a jet 'redirection' effect (see beio.v;. 
conclusions were reached; 


number of Impinged 
required to remove 
The theoretical studies 
to model it explicitly 
The f o 1 1 ow i na 
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o By simplifying the floor flow, tangential blowing ;s likely 
to have eased the task of modeling Impingement cases. 

o The use of a 'kinked' theoretical jet plume, with line elements 
along the tunnel floor, may cause spurious Interference effects, 
particularly if vortices are used in the flow model. 

o The Introduction of a theoretical 'plume redirection term 

(Figure 5.11), which restores the plume from it's kinked forn to 
a free air trajectory, can reduce the apparent interference by 
a factor of three or four. 

The last finding is considered very significant. 


V 


28 



ReFERENCES 


1 . 


2 . 


3. 


I4. 


ream R. L. and Weston. K, P.. ’•Velocitv Field of a Round Jet in a 
CrossVlow for Various Jet Injection Angles and Velocity Ratios. NA 
Technical Paper 150b, Oct,, 1979. 


Heyson H. H. , "General Theory of Wall Interference for Static 
Tests in Closed Rect^nqular Test Sections and in Ground E. ect 
Technical Report. NASA TR R-36^, Sept. 1971. 


St ah i 1 i t V 
" NASA 


Williams J. and Wood, M. N., "Aerodynamic Interference Effects with 
Jet-Lift’v/STOL Aircraft under Static and Forv^ard-Speed Conditions. 
Zeitschrift fur F 1 ugw i ssenschaf t en , Vol 15. Sec 7, pp. 237-256. 

July 1967. 


Weston, R. P., "Vorticity in the Jet in a Crossflcnv, 
to the AlAA South Eastern Regional Student Conference; 
also NASA TP 1007 "Induced Velocity Field of a Jet m 
1978, by Fearn and Weston). 


Paoer G9 , presented 
April 1973. (See 
a Crossflow," May 


29 



• suonejnb I jijoo (puojjuaAuoo jo^ uiejbojd 

^.sAieue ajnjeubiS.ajnssajd-t|P^ ‘l l-'^d 1 said.r^uiJd bujje-iado r2 ^^nbtj 


NftON^ woyd 

3iVmtfA3 


3NIWy3l3Q 


iiidino 


J WnON>t W0b3 

3ivrnvA3 


S3ivn3ais NO 
t) 30 133333 3Hi 
a03 S13303 33N3m3NI 


S3«niVNU»S 

nvft3Qis 


3niwyiJ.3u 


snvn3QJS 
NO 3 30 
133333 3A0W3« 


S3bniVND»S 

yooi3 ‘i 30oy 



1 ■— 1 


S3anss3yd i3NNni 03ynsv3w 


yoo33 ^ 300y nc 
J 30 133333 3H1 
y03 S13303 33N3m3NI 



S33Q0W n013 ltf3ll3y03Hl 


V 









•sujpjbojd I jjPd pue juasajd atjj ua.^Miaq a^ejja^ui z'Z 


(l*Z 3bflDU 
33S ) 
WV'dDOad 
i idvd 


uvdoodd 

lN3S3dd 












































BLOCKAGE 

RATIO 

Uc/Uo 

.20-r 


.16 + 



1" JET 




JET 




Figure 3-5 Supply pi 


TUNNEL q_ 


— 30° 



36 










Figure 3-9 Nominal velocity ratio as a function of plenum pressure ratio. 

UQ 






0 1 2 3 ^ 

q 

NOM 


Figure 3-10 Flow coefficients for 
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Figure 3. 15 Wal 1 pressure 


orifice locations for jet - i n-crossf 1 ow experiments 













b) Downwash distribution 


Figure ^.5 E><it velocity contours for the 1-Inch jet 
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3” Jet . pipe prof i le 


Figure A. 7 Effect of jet si 
a) R » 2 
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Figure ^.7 (Continued) 




igure ^.7 (Continued) 
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1, 8 Vortex location as a function of jet velocity ratio 
a) Effect of jet diameter and profile shape 
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Figure ^.11 


Vortex-lift correlation for 1-inch jets 
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Figure 4.12 Jet pipe datum runs at RnOM = ^ 
(a) 3-inch jet pipe 
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Figure ^.12 (Continued) Jet pipe datum runs at R^OM = ^ 
(b) 3-Inch jet pipe with 'ground' BLC 
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Figure 4.12 (Concluded) Jet pipe datun runs at R = 0 
(c) 1-inch jet pipe 
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Figure 4.15 Wall pressures for 3-inch jet 'square' profile 
(a) ^nOM ~ 2 
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Figure ^ 4 . 1 5 (Continued) Wall pressures for 3-inch jet 'square' profi 
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Figure ^.15 (Continued) Wail pressures for 3-inch jet 'square' profile 
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Figure 4.16 Wall pressures for 3-inch jet Pipe prof i 
R * 6 and R * 8, with 'ground' BLC applied 
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Figure ^.19 Wall pressures for vertical and inclined jets at R = 8 
(a ) Rails ^ and 3 
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Figure . 1 9 (Concluded) Wall pressures for vertical S inclined jets at R = S 
(b) Rails 2 and 1 


6 




SOURCE/ 

DOUBLET 

PAIR 


VORTEX 

PAIR 














30-i 


Z/d 


20H 




04 — 

0 


— fearn curve 
FOR JET CENTER 

— WILLIAMS 
6 WOOD 



Figure 5-1 


(Concluded) Flo. ^del geonetrv 
(c) Selection cF trc-ijectorjes 
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Figure 5.3 Wa!i pressures for two available flow models 
(a) The "Fearn" , vortex model. 






Figure 53 (Concluded) Wall pressures for two available flow nradels 
D • 3" . R * 2 . 

(b) The "Heyson", lifting-doublet model. 
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5.U Comparison between jet-ln-crossf low models, 
(a) Present, 'VSO' flow model. 
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5.U (Concluded) Comparison between jet-in-crossflow models 

(b) Heyson's flow model. 
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WAKE TRAVERSE DATA FOR ONE- INCH JET 
AT -30°/ 0° AND +30° TO THE VERTICAL 
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Figure A1 (concluded) Effect of jet inclination on total 
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F igure A2 (cont i nued) Effect of jet Inclination on vertical velocities: 
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Figure A3 Effect of jet inclination on lateral velocities: R~2 
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Figure A4 Effect of jet inclination on streamlines: R»2 


1" Jet , 30- forward 


I" jet, 302 aft 










-.12S 
-.10 
-.08 
-.06 
-. 0^4 
-.02 
.00 
.02 
.Oi* 
.06 
.08 
.10 
. 12 S 


1“ Jet, 30- forward 


Figure A4- (concluded) 


■''J I I 
X / 



1 '* Ver t i ca I jet 


I I 

\ \ 

\ \ 

\ \ 

V 


\ 


-.u 
-.10 
-.08 ^ 
-.06 
-.04 


I ^ 

^ II 
/! 

^ // 

/ 

y 



-.02 / 
.00 / 
.02 


.04 

.06 

.08 


. 10 
.14 


1" j' 


Effect of jet inclination on streamlines; R*8 




APPENDIX B 


WAKE TRAVERSE DATA FOR ONE-INCH/ 
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Figure B1 Effect of jet profile and size on total pressures: R=? 
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Figure 81 (continued) F.ffect of jet profile and size on total pressures: R-6 
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Fiaure B2 (continued) Effect of jet profile and si?e on vertical velocities: P=3 
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Figure B2 (continued) Effect of jet profile and si^e on vertical velocities; R=^ 
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Figure B2 (continued) Effect of jet profile and size on vertical velocities: R = 6 
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Figure 03 (concluded) Eff-ict of jet profile and size on lateral velocities: R = 8 
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Figure B4 (continued) Effect of jet profile and size on streanhnes: R 4 
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Figure e<l (continued) Effect of jet profile and size on streamlines: R-6 
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Figure B4 (concluded) Effect of jet profile and size on 
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APPENDIX C 


SOME comments upon HEYSON'S METHOD 


Heyson's model differs from the VSO model not only qeometrlcally 
(see Figures 5.1(a) and Cl) but also because it relies entirely upon 
doublets. A further significant factor is that, as a result of being force- 
based, the form of dependence upon R may differ, 

Fvoff'^ties of line doublets 

Heyson's model involves only line doublets. However, these can have 
source-like properties which permit for example, a non-zero, wake-b 1 ockage- 
1 i ke asymptote far downstream. The simplest example of this concerns the drag- 
related line doublet which streams along the tunnel wall from the impinge- 
ment point. The upst ream-d i rected doublet vector aligned with this line 
element becomes in effect, a long thin ring-vortex tube which pumps fluid 
along it's center in the upstream direction. It may be shown mathemat i ca 11 y 
that the induced flow field is identical to that for a source situated at 
the upstream end plus a sink far downstream. The volume flow Q, ftVsec is 
numerically equal to the line doublet Strength u ft^/sec/ft. The inclined 
doublet line from the jet exit may be regarded similarly, giving a source at 
the upstream end and an equal sink at the impingement point each with a 
strength proportional to the doublet strength component parallel to the 
line between them. However, these contributions create the equivalent of 
positive solid blockage (rather than wake blockage) since the net source 
strength within the test section is zero, A vertical doublet component, 
cor respond i ng to lift, has no wake-like contribution. it has a negative 
solid blockage effect when considered as above. 

The behavior for ang 1 e-of-at tack interference is the converse of that 
Just described: lift-related doublets generate a positive A.t asymptote 
downstream and drag-related doublets give an antisymmetric pattern of 
urwash and dcwnwash. 

Though the individual elements have the properties just ment i oned : i . e . 
appropriate asymptotic behavior and the ability to provide 'peaky' behav<or 
near the model, the geometric restrictions of Heyson's model are quite 
severe - particularly with regard to the limited range of jet skew angle, 
which is found 'n practical cases. (For 2 ^ R 8, 58,3'’ ^ f A8.6 where 
Y is from the vertical. 


DependeKOu upon Jet-velooi tu 
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The form adopted by Hevson is a genera 1 i zat i on of the familiar 
5 , S/C ‘Cl equation to include cross effects. Thus 



: 

f, 


Uo “ 

"11 

3.1 

;,w 

A 

f 


u« 

■ 13 

1 

3.3 


(Cl) 


156 



where 5ii and 5-j-i are the Influence factors for drag on u and lift on w and 
fi and fi are related to drag and lift respectively. These take the form of 
the product of mode 1 -to-tunne 1 area ratio and a characteristic velocty. Thus 

° (C2) 


T X S 


T 


Uq and Wq are defined via 


0 - and L - - cA„ 


where the product pA„ V„ represents the jet mass flow In the present case. So 

D - - m Uq and L ■ * m Wq 
giving Uq * ■ D/m and Wq “ * L/m 

If A,^ is taken here as jet exit area, then Vr - Vj and A^^ - m/. Vj 
Substituting into Equations (C2) from (Ci*) and (C3) 9' 'as 

m/pVj 



-D 

(TH 


and 


-L 


m Us» 


(C3) 

(CM 




and 


(C5) 


Li ft 

L is the net lift on the model, which we shall assume equals jet gross 
thrust for the present vertical jet. Thus 


Vj 




(C6I 


giving f rom (C5 ) 

^3 ‘ pVj i-uC ’ C 

This is the form needed if Heyson’s program is used directly. For the 
calculation of wall pressures, however, a special version of the present 
program is used which reproduces the Heyson model. Thss requires input of 
doublet strength, given by 


'j Vj' 


2- o 

1 


2 " 


from (C3) 




(C7 
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Drag 


For the most meaningful application of Heyson's method to the present 
case it would be preferable to employ balance data, In it's absence, jet 
drag will be estimated from the traverse data ar X/0 - 6, since experimental 
vortex strengths and core sizes are available, tor a solid-body vortex core, 
which is a reasonable approximation to the present case: 

where b is the vortex spacing and a is the core radius. It was noted from 

experiments described in Section h.} that the vortex cores touch, implying 
that b ■ 2a. For this case, the argument of the log in Equation (C8) is 
unity cirid the first term becomes zero. Thus 



(C9) 


Substituting Equation (C9) into Equation (C5) 


Now 
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Th^ R- cubed power in Eq-.ation CCIl) is unexpectedly high However 
it w::Te. ,han^his (,.oug. .tU, ,rea.e. ,.a. «-s,.a.ad) 
because the relationship for b takes the emp,r,cal forr, 


b = (a + 6R) 0 


(C 12) 


, rri Q - n 1 for the 1-inch let dat at X/D - 6 in Figure h.Sia}. 

where cj « 3 aod 6 -3 into (CM) that the resulting values of 

It is found upon suostitutmg (C12) into (CIM tnat tne '^e^u y 

u are far too small; the implied cylinder d.ameter values are 

. . , a- fhe let The assumption was therefore made that b » U- 

magnitude less than the jet. ne a^bumpi of the cubic 

thi.; nave a match to the jet diameter for R - 7.33. Because or tne cuo c 
depenLncy (Equation (Cll ), the effective diameter decreases very rap.div 
wUh R ^n making the substitution b - D and clearing, we obtain 


_L 

32 


(C13) 


and 


1 


U 


1 R 

TT 7T 


Unless stated otherwise, Equation (C13) O'" i^of this 

doubLt strength for the Heyson examples quoted m the mam. text of th.s 

report . 


Jet ske-j angle \ 

The sweep-back of the jet doublet line, from the vertical 
to the velocity ratio R by the equation 


is related 


X 


j * tan ' 


R 


(CIS) 
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PROGRAM description 


Cav^hilitics 

This program computes the effect of a single jet of an arbitrary 
diameter at an arbitrary location with respect to the global, tunnel 
coordinates. Using the model described in the ovain text and using an 
itnaging technique, the additional jet effects are calculated at the spec- 
ified wall points and at the tunnel centerline. 

A limited option is provided via the flag ICAL (see Input Description 
below) to study the effects of individual components of the theoretical 
nodel. if so desired. An option flag JMOD also lets this program be run 
to study the effects of jet only when it is not desired to run the Part I 
program for further analysis using the modified wall signatures. 

As of now the theoretical jet model is restricted to jets issuing at 
90“ to the mainstream. To consider the effects of jets at other angles, 
the input constants as well as the program coding defining the geometry and 
strengths of singularities will have to be appropriately modified. 

I' the Part-I Prcji^-V-> 

To ot.ain the total interference of a model with a lifting jet, the 
Part 1! program inist be run first as a pre-processor. This program creates 
two mass 'torage files via FORTRAN UNITS tO and 11, which subsequently 
beco™ part of the input files for the Part 1 program, with the same 
FORTRAN UNIT-numbers. Additionally, the Part I program must be "signalled' 
to expect the pre-processed jet-effect output, This is done by assigning 
a non-tero value to a variable JETEFCT, which is the last variable of 
Input Card Number-2 in the Part I program. 



Inj'ut Desarirtion 


1 I I CAL [ JMOD 
ICAL: 


For/T« t 1615 

A calculation index (0 < ICAL ai 

0 - No calculations are to be performed. 

Print the model geometry and singularity 
strengths . 

1 - Calculate jet effects using vortex pair 

singular! ty only. 

2 - Calculate jet effect using source singularity 

on 1 y 

3 - Calculate Jet effect using doublet s i ngu 1 c. r i t y 

only. 

k “ Use all three singularities (Recommended value^ 


jmoD: a non^zero value implies that the effects of jet 

should be taken out of existing wall pressure 
signatures. The wall pressure signatures to b-.'^ 
modified should be available in UNIT •* 7. The 
modified signatures will be written to UNIT - 10, 
and the interference effect due to jet alone will 
be written to UNIT - 11. 


If JMOD » 0, no input Is sought from UNIT - 7. The 
wall pressure signatures due to the jet alone are 
written to UNIT - 10, and the interference velocities 
due to the jet alone are written to UNIT - 11. 
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Format {20A^) 

A title with no more than 80 characters. 

Format (8F10.6) 

Tunnel width 

Tunnel height in same units as that of 8. 
Format ( 8 IO. 6 ) 


AZ 

BZ 

CZ 



AY 

BY 


Non-dimensional constants For definition 
Z-coordinate of vortex curve (see Eg. 5.1K See 
sample input for recommended values. 

Format (8F10.6) 

Non-dimensional constants for the de‘’'nition of 
Y-coordInate of vortex curve. {See Eg. 5.3) See 
sample input for recommended values. 


6 


AG 


GB 


CG 


Format 


(8F10.6) 


AG Non-dimensional constants for the definition of 

BG circulation strengths. (See Eq. 5.^) See sample 

CG input for recommended values. 


7 


ACL 


BCL 


CCL 


Fornat 


(8F10.6) 


ACL 

BCL 

CCL 


Non-dimensional constants for the definition of 
Z-coordinate of jet centerline. (See Eq. 5.2). 
See sample imput for recommended values. 
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Forrnat (16I5) 


NSPV 


NCAL 


NWALL 


LAYER 


Number of singularity points In the X-dlrection 
on the vortex line. (Note: These points will be 

connected by straight line segments to produce 
the required links. The same number of points 
are used fo* sources and doublets) 


NCAL: 


NWALL: 

LAYER: 


Number of calculation points in the X-direction 
for both wall pressures and the centerline inter- 
ference velocity points. 

Number of rails on tunnel surfaces. 

Number of image layers to be used. (Recommended 
value: 5) 


XNG, 

XNG 2 

• • • 

XNGnspv 

XNG. : 


X 

-coord i nate 


I 


with respect 


Format (8F10.6) 

of i-th singularity point normalized 
to tunnel breadth, B 


XC, 

XC 2 

... 

XC 

^NCAL 


Format (8F10.6) 


XC. : 

I 


X-coordinate of I-th calculation point normalized 
with respect to tunnel breadth, B. 


YWALL ^ 

YWALL 2 


"“''‘■'■NWAU 


Fermat (8F10.6) 


YWALL. : 

I 


Y-coord I na tc of i-th rail on tunnel surface 
normalized with respect to tunnel breadth, B. 






12 


ZWALL ^ 

ZWALL^ 

• • • 

^^ALL^yALL 


Format (8F10.6) 


ZWALL. : 
» 


Z-coord I date of i-th rail on tunnel surface norm- 
alized with respect to tunnel height, H, 


13 


IFLOR 

IROOF 

1WAL1 

IWAL2 


Format — (l6l5) 


IFLOR: 
IROOF. 
IWALl : 
IWAL2: 


The rail number for floor signature. 

The rail number for roof signature 
The rail number for the signature on sidewall 1. 
The rail number for the signature on sidewall 2. 
(Note; These should be compatible with the values 
defined in Card No. 10 of the input for 
Part I program) 


lA 


ITEST 


IRON 


(POINT 


Format ( l6 1 5) 


ITEST: Test Number 

IRUN: Run Number 

I POINT: Point Number 


m 

D 

YJET 

ZJET 


Format 


(8F10.6) 


R: 

D: 

YJET: 

ZJET: 


Velocity ratio, 

Jet diameter, normalized with respect to tunnel 
height, H. 

Y-coordInate of jet origin normalized with respect 
to tunnel breadth, B. 

Z-coordinate of jet origin normalized with respect 
to the tunnel height, H. 
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Mass Storage Requirenerits 

In addition to the standard Input and output files (Uh41TS 5 and 6), 
the coding also requires that three more mass storane files be pre-assigned 
to UNITS 7, 10 and 11. The purpose and format of these data files are as 
fol lows : 

UNlT-7 Input file. This file should have the measured 

wall pressure data for the tests done with the 
model and the jet. The format of the data is the 
sane as the corresponding wall pressure data 
file used by the PART I program. This file need 
not be assigned if the jet alone option is used 
(i.e., JMOD » 0) in running this program. 

UNIT-10 Output file. This file will contain the wall 

pressure signatures as modified by the presence 
of the Jet. This will be in the same format as 
UNlT-7 and can be used as it is for the auxiliary 
input file for Part I program. If the option 
JMOD ■ 0 is employed, the file will contain the 
wall pressures due to the jet alone. 

UNIT-11 Output file. This will contain the interference 

velocity at the tunnel center 1 i ne due to the jet. 
The Part I program will add these values to the 
model interference to obtain the total interference 
due to the model with jet. This file should be 
assigned to UNIT-11 again while running the Part I 
program for further analysis. 
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Output description 

The output produced by this program consists of the following 
sections: (1) Tunnel geometry and the values of constants used for 
defining the model as specified in the input to the program. (2) For a 
given jet diameter and R-value, the complete definition of the model. 

(3) The velocities induced by the jet at the pressure rails on tunnel 
surface and the new values of wall pressures as modified by the jet. 

(A) The interference velocity due to the jet alone. 

A listing of the program, a sample Input and output are given below. 
The following list of subroutine^ are the same as in Part-1 program and 
they have to be additionally linked with the present program: INFLU, 

CRDTEG, LNVXGN, LNVXEQ, LNSCGN, LNSCEQ, LNOBGN, LNDBEQ, CPREAD, 
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non 


PROGRAM OETFFFCT 


C 

C 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

r 

C 


FILE... JETEFFCT . 'OR ; I 


Ca1cualt*3 effect of a Jet on tunnel w« ’ 1 » and modifies the wall 
,ccorV).ngly. d.,ter™ ne, II,,- l„t,rf.,.nc 

velocity conponenls at the tunnel tentet^l'nr cue tie Jet. 


IhPUT flies ! 

U»IT_5 

OUTPUT Files: 

UN’T_6 

'IN n_,i0 

■JN : T_i I 


Flow parameters and oti er key values. 
Meisured Wa 11 _p r e:. sur e dsts. 

Print O' ' 

M.^difled W »1 l_orensure datt. 

Int. crferenre vslcn it’s me to the let 


r 

C 


DIMENSION XPV( 30 ) , YPV( 30) , 7PV1 30) .XPR I 30 ,YPR 

I XPSf :u ) . YPS( 30 ) .7PS< 30) . kC(30 . VC 

I VMS < 3 H ) . GAM! 30 ) , VMU ( 30 > . 11301 V { !0 

I liWL I 30 > , I'Kl 1 30 ) . VVAL L ( 10 ) . 7VAI I, ( 10 ) 

COMMON /IMAC/ 8 ,M. layer .mini T . IBU'. 

CCHMON /DIPC/ DCXO. OCYO, [’CZO 

CCMMON /LOAU/ AL f U . POU . 'JOU . C MUU , C L U , C 3 U . ( MU 


3f > ,7PP ( 30 > .XNGl 30) , 
3£ > . 7f ' 30 ) . 

■ . V ; 30 / . i:f' 1 30 . 3 ) . 
n VP { t : . T 1 TL ( 7 j1 > 


PI * 3 . 141597^54 
IBUG = 0 


REA0_.IN GEOMETRICAL FAPAMETERS FOR THE JF^_M0DFt 


r 

t, 

C 


RFAD < 5.5 501 ICAL.JMOt’ 
READ (5.570) TITL 
REAC (5.500) H.H 
WP 11 E ( 6 . 600 ) e . H 
'J? n E ( G . 60 1 ) TITL 
PfAD (5.500) AZ.BZ.cz 
WR I I F. <6.61?) 


WP I 1 E (6.630) 
RLAD (5.5ft0) 
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